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Abstract 

We describe a generator for the process e + e~ — ► W + W~ including all 
one-loop and leading log photonic two-loop contributions. It includes po- 
larization of the beam and W bosons, and the possibility to study the effect 
of anomalous couplings. 



1 Introduction 



The computation of the one loop W pair production cross section in e + e~ colli- 
sions has been performed a long time ago [I], [^, 0]. The effects of hard radiation 
have been added to the last two calculations f|, [5]. These results agree at the or- 
der of 10~ 3 or better ||, except at high energies in the forward direction, which is 
probably attributable to numerical instabilities. With the advent of LEP II there 
is a need for an event generator for W pair production at this level.Q We present 
here such a generator based on the computations of Ref. 0. It includes the full 
one-loop matrix element, leading logarithm two- loop initial state bremsstrahlung, 
polarization effects and anomalous couplings. The main limitation is that we as- 
sume the W bosons to be stable particles. Work to remove this limitation is 
in progress. We nevertheless think that the present version is useful in several 
respects: it includes the interesting part of the full process e + e~ — > W + W~ — > 
four fermions and thus allows to study what we may learn from VT-pair produc- 
tion under idealized circumstances. The missing part only involves well known 
physics, namely, the decays via the well established charged current interaction. 
The decays thus just serve as polarization analyzers for the W s. Since a full 0(a) 
calculation of the observed process e + e~ — > four fermions is a major enterprise, 
we think it will be useful for tests of parts of more complete calculations. 

We first recapitulate the construction of an event generator, then discuss the 
polarization, initial state collinear bremsstrahlung and anomalous couplings. In 
the appendices technical details about the installation and sample input and 
output are given. 

2 Event generator 

The event generator presented here is based on the calculations 0, [7|. These, 
however, have all been presented as total and differential cross sections, or cross 
sections with certain cuts in the photon energy and the angle between the photon 
momentum and the beam axis. In the present work we present an event gener- 
ator, i.e., a program which will give configurations W + W~ (7) with probability 
proportional to the contribution to the total cross section. This allows experi- 
mental studies to be done including arbitrary cuts and detector capabilities. In 
this section we describe how the conversion was performed. 

The last step to an event generator is simple: if the cross section is written 
as an integral over a positive function f(x), with x a point in a hypercube < 
Xk < 1, and / max the maximum value of the function in this hypercube, event 

1 After writing this article it has come to our attention that a similar effort has been reported 
before f§. 
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generators are given by the representation of the cross section 



o = jdxm^Wm^^Y.v^-^), (i) 

i=l i=l J max 

where the £j are random points in the hypercube and the random numbers 
between and 1; N is the number of function evaluations used. The first sum 
corresponds to a weighted event generator, the second one to a weight-one event 
generator — all events correspond to the same probability. Both options are incor- 
porated in EEWW, the first using the adaptive Monte Carlo routines vegas ||, the 
second one using non-adaptive algorithms simplemc and axmc. Experimental cuts 
and efficiencies can now easily be incorporated in the sums as ^-functions, pos- 
sibly involving more random numbers. For parton-level calculations a weighted 
Monte Carlo is generally more efficient, as the main cost is the evaluation of the 
matrix element. When a detector simulation is included, however, this part will 
take much more time, so a weight-one generator is generally used. The efficiency 
of the weight-one generator is given by 

t-fv&vjihm-r,), (2) 

J J max i = i J max 

i.e., the average number of events generated per function evaluation. 

The maximum value of f(x) may not be known beforehand. We obtain an 
estimate by sampling the function at 10 1+n / 6 points for an n-dimensional func- 
tion, and taking 20% more than the largest value found. This is adequate for a 
reasonably flat function. If during the generation a function value is found which 
is larger than the assumed maximum value, Eq. (P is adjusted as 

efip. - ,,) _ int M) + a (Ij#L - lnt ZM - ,,) , p) 

/max /max /max /max 

where int(a;) denotes the largest integer smaller than x. This means that the 
event is accepted as many times as is necessary to obtain the correct integral. As 
a result, not all events are independent, and statistical fluctuations may be larger 
than expected. Note that / max is not readjusted to reflect the new maximum 
found. 

For our purpose we need to write the one loop cross section as an integral 
over a positive function on a hypercube. The hard radiation is already in this 
form: 

erf = — / d(f) W / dcos6 w / gLEv / d0 7 / d cos 9~J\M? | , (4) 

2s JO J -l Jk E e JO J -I 

with J the Jacobian defined in Ref. and E e the beam energy E e = \/s/2. 
kg is the minimum fraction of the beam energy the photon is allowed to have. 
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The mappings from the X{ to the variables in Eq. (B) are also described in this 
reference; these map away the infrared 1/E 7 , the collinear l/(E e ±p e cos # 7 ), and 
the t-channel 1/t peaks of the matrix element squared. 

The Born, virtual and soft contributions to the total cross section have the 
same form without photon integrals, so the integral can be taken the same with 
a flat behaviour in these dimensions. These two integrals are combined by intro- 
ducing an extra (sixth) integration variable 

o = a + ^ V+S + *? = /&* (* + a^) + ^^a*} , (5) 

jo L a v J 1 — a ) 

with < a < 1 a tunable parameter. The distinction between hard and soft 
radiation is arbitrary; the only demand is that the cutoff k Q is much smaller 
than unity to validate the eikonal approximation used in the soft bremsstrahlung 
integrals. We use this freedom to make the first term as small as possible without 
making it negative; a suitable choice is 

7T 7T 2 3 , 2 1 \ 

h ^ S ~ A ~2^~Y + 1 ~4 l0gK/g) , R , 
fc ° = ^^ eXP \o g (K/s) - 1 • (6) 

The parameters A = 130 GeV, B = 35 GeV in the a scheme (B = 9 GeV 
in the scheme) and a = -y/i/1000 GeV below ^Js = 350 GeV, a = 0.35 
above have been found by trial and error. At LEP II energies this gives (<To + 
a i +S )/ cr f ~ V^- As the virtual matrix element takes much longer to evaluate 
than the bremsstrahlung amplitude it is advantageous to keep this ratio as small 
as possible. Using these values we obtain a reasonably flat function f(x). The 
efficiency, without any cuts, is about 40%. This translates to about 5 events per 
second on a workstation. 

Note that the resulting soft photon spectrum, which seems to continue down 
to below 1 MeV, can not be trusted below a few times the W width because of 
the narrow width approximation. 



3 Polarization 

Building on the calculations presented in Ref. |7j the event generator has many 
possibilities for studying the effect of different polarizations of the initial and final 
state. 

The state of the beam is characterized by the electron and positron density 
matrices in the helicity frame 

_( I + Pl Pre-* \ (7 , 
p - ~ 2 ~\P T e^ 1-P L ) [{) 

i + ff-P' _( l + PL P' T ^' \ (R) 
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where P = (Pt cos 0, Pt sin 0, Pl), Pl is the longitudinal polarization and Pt the 
magnitude of the transverse polarization of the electron; P' L , P' T of the positron. 
The polarization vector of the electron (positron) is pointed in a direction ((f) 1 ) 
upwards (downwards) from the outward direction. For the natural polarization 
in a storage ring = 0' = 90°. Special cases are no polarization (Pl = Pt = 0) 
and longitudinal polarization (P T = P T = 0, P L = ±1, P' L — for right- and 
left-handed polarization); these can be chosen separately. In the limit that the 
electron is taken massless the amplitudes for like-handed electrons and positrons 
(Pl = P' L ) are zero. The matrix element is then given by 

iw = ^[(i-^D(i^ + i 2 +i^-i 2 ) + (^-^)(i^+i 2 -i^-r) 

+ (2P T P T ) (cos(0 -(f)' - 2<f) W ) Re(M + M*_) + sin(0 - 0' - 2cf) W ) Im(M + M*_)J§) 

A detailed discussion of the effect of the choices is given in Ref . .[] 

The W bosons can have three polarization states: 2 transverse and 1 lon- 
gitudinal. One can either choose to average over these in the matrix element 
('unpolarized' Ws), or generate them separately ('polarized'). In the latter case 
the polarization vector is made available, so that the subsequent decay of the W 
can take this information into account. The other two possibilities are to generate 
only longitudinal or transverse bosons. 

The generation of different polarization states is implemented by converting 
the sum over final states into an integral as 

1^ o(\ + ,\ ) = dx}_^ — a(A+,AJ, (10) 

A±=-1,0,+1 i=0 Ci +! C% 

with = int(z/3) — 1, A~ = mod(z, 3) — 1. The cutoff values q are chosen 
according to the lowest order matrix element squared at the same W~ angle: 

|M (A+ Am 2 

Cl+l Q "E? = o|M (A;,A7)r W 

The integration variable to make the discrete choice ( |T0"D is taken to be the 
rescaled version of the one used for the choice between the hard and other cor- 
rections (0). 



4 Initial state collinear bremsstrahlung 

The largest corrections in the one-loop calculation are caused by initial state 
collinear bremsstrahlung. It is therefore natural to seek a way to include these 

2 Note that in table V of this reference Pt = Pt = 1 is taken, like in the definition (21) but 
unlike the suggestion in the text, and a factor two has been included in the definition of the 
asymmetry. The typo corrections coscj} — > <fi in (21) and A0 = </> — </>' — 2(f>w after (17) are 
obvious. 
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to higher order. The leading log terms are well known [BJ. They were included 
in Ref. || by replacing the lowest order term cr (s) by 

a im (s)= / dkp ini (k)<Ja(s(l-k)), (12) 



o 

with the function p ini given by 

p- mi (k) = (3k^(l + S^ S + 5^ +s )+Sf + 5 2 H , (13) 

with (3 = 2^(L — 1) and L = log(s/m 2 ) the large collinear logarithm. The terms 
S denote the infrared finite parts of the leading one- and two-loop corrections. 
Explicit expressions are given in Refs |J |j, except that we use Sf = -((L — 
l)(fc — 2) + k), which differs by the last (non-leading) term.Q However, one 
has now included the one-loop leading corrections twice; this is corrected for by 
subtracting the unexponentiated 0{ot) contribution from this formula: 



^.double 



logfc + <^+ S ) <ro(s) + / fc femaX dk(£ + Sf)a (s(l - k)) , (14) 



with k Q an arbitrary (small) cutoff. 

This procedure, however, is not suited for an event generator. First, a separate 
structure function has to be assigned to each incoming particle [0, [II], [L2|, |i~3H ; 



these are identical to p ini with j3 — > /3/2; for simplicity we will keep the single 
integral notation. A more subtle point is that the structure function ([13|) has been 
computed by solving QED evolution equations after integrating out the angles of 
the extra photon over all phase space. This inclusive approach is not appropriate 
for an event generator. We solve this by dividing phase space into two regions: 
inside a cone of opening angle 9 C around the beam direction one declares all 
photons collinear (and makes an inclusive measurement), outside it one measures 
the angles (and measures exclusively). Only in the phase space inside this cone 
(9 C is assumed to be much larger than m 2 e /s) we use the exponentiated two- loop 
structure function. Outside this cone the Oipt) result is used, with the possibility 
of extra collinear photons. The two limiting cases are 8 C = 0, which gives a strict 
0(oc) behaviour, and 9 C = n, which reproduces the totally integrated cross section 
computed before.^ Integrating only over the angles inside the cone amounts to 
replacing s by a scale p 2 = s(l — cos(9 c )/2 in the collinear large logarithms: 
L = log(/i 2 / 'ml) . A calculation of the px dependence of the structure functions 
has not been done yet and is beyond the scope of this article. 

A technical problem is that adding and subtracting parts does not give the 
positive definite integrand needed for an event generator. To obtain this we 

3 Numerically, the difference is less than 0.1%. 
The region ir/2 < 9 C < it is redundant from an experimentalists' point of view, but the 
collinear terms have to be integrated over all angles to obtain the full result. 
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rewrite the total cross section as 



a 



a 



a 



v+s 



a 



double 



dkp ini (k) cr (s(l - k)) + 



a 



V+S/ 



((3\ogk + 5^ +s )a (s) 



+ 



Mi 



dktf- + 5 H )a (s(l-k)) 



(15) 



The terms in square brackets are of order a and do not contain any large collinear 
logarithms; we can therefore add them to the first integral without introducing 
terms of order a 2 L: 



a 



dkp ini (k) (a 



Z.V+S 



a 



(16) 



with a\ the differences in square brackets in Eq. (fT^). The virtual and soft one 
is easily computed, to subtract the leading log from the hard radiation we have 
to reintroduce the cos 9 integral as (see, e.g., WM) 



dk 



ko 



(1 
k 



a 



7T Jk E e 




cos 8 r 



dcosO 



Pe 



E e + p e cos 9 



Pe 



dcos9- 
cos e c E e — p e cos 9 



(17) 

This can be subtracted from the hard radiative integral over a 3-particle phase 
space if one adds a <fi integral (which is trivial for Pt = 0, otherwise the same 
0-dependence is taken as the lowest order). A similar approach is followed for the 
double pole terms; here the upper integration boundary is of order m\/s(l— cos 9 C ) 
which can be neglected when 9 C m 2 e /s. The final result is just the collinear 
limit used already in Refs [|1^, |!| , but now applied for all angles up to 9 C in both 
the forward and backward direction. 

Unfortunately there is no reason for the integrand of of to be positive definite; 
indeed, in the collinear limit it is zero. Recalling that within the collinear cone we 
do not make a distinction between different collinear photons we add a fraction 
of the Born cross section to restore positivity. We take 

dk( + )dcos9( ^— 

71 Jk V-l JcosOj v k 

-j^jdob^l - k)s) da (s 
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P±-Pl 



+ (l-k) 



n 



=v+s 



a 



a 



X 



v+s 



d cos 9 d cos 9 

- 1) log A; 




max ~ ( 2 2fc max 



7T 



)L H — 2 + 2A; max 



h 

_ run 



where p± is the positron (electron) momentum in the first (second) cos ^-integral. 
Note that /c max depends on the angle between the W~ and the incoming electron. 
As the expression for the hard bremsstrahlung is dominated by the added Born 
term we generate an event with only collinear photons in this region. 



5 Anomalous couplings of the W boson 



One of the main purposes of LEP II being the study of the WW^f and WW Z 
interactions we allow for the inclusion of non-standard couplings at these ver- 



tices [16, 17]. We parametrize these by the CP-invariant effective interaction 



Lagrangean for the WWV interaction (V = 7 or Z) 

4m = ic v {v»(w+w-»-w-w + ») + v^wtw-} (20) 

A£ e y ff = iC v S v {v v (W+ v W-» - W-W +>1 ) + V^W+W-} + iC v 5k v V» v W+W; 

+ ic v ^-v^w + ^w^ + Cv^^v"") {wfd v w; + w-d u w+hi) 



Electromagnetic gauge invariance forces 5 1 = 0. This effective Lagrangean is 
related to the other widely used form |18|] by g\ = 1 + 5y, fty = 1 + Sy + 5k v 



and fcf = £v s / m w- Whereas CP- violating terms are known to be small, we see 
no reason to neglect the parity- violating form factors £y. 

The inclusion of these anomalous couplings in a one-loop expression causes 
problems, as the resulting theory is non-renormalizable. We assume that the 
anomalous couplings are small perturbations of the standard model, and that the 
leading effects enter only in the Born term. They are therefore not included in the 
0(a) corrections, in particular not in the hard (non-collinear) bremsstrahlung. 
The usage of a different expression for the Born term in the soft bremsstrahlung 
will upset the validity of Eq. (P) which assures the positivity of the integrand. 
This only causes negative events when the anomalous couplings decrease the Born 
matrix element squared, which, due to the violation of the gauge cancellations, 
happens only in the backward direction (cos 6w ~ —1)- We therefore increase 
the cutoff by a factor 1/xf, where 25 is the random variable which is mapped 
to cosdw- In case of problems one can rescale ko by an arbitrary factor. At 
high energies, for large anomalous couplings and with extra initial state collinear 
radiation it may not be possible to avoid negative points. 



6 Anomalous couplings and loop effects 

We have compared the effects of anomalous couplings of the WW Z vertex to the 
one-loop corrections in the total cross section and the angle of the W~ boson. 
One should note that the effects of the one-loop corrections is to entirely different 
amplitudes than those affected by anomalous couplings. The effective one-loop 
contributions to the anomalous couplings defined in Eq. ( pip are of order a/71", 
i.e., a few times 10 -3 . This is completely negligible at the precision expected 
at LEP II. The large one-loop contributions therefore originate in different form 
factors, like the ones associated with the t-channel graphs. These also affect 
simple distributions like da/d cos 9 w- 
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The values used for the parameters were mz = 91.176 GeV, mw = 80.152 GeV, 
m t = 174 GeV, m# = 100 GeV. We assumed unpolarized beams of 87 GeV each 
and used the renormalization scheme. The results are shown in Fig. [1]. One 
sees that the main effect of the radiative corrections is described fairly well by a 
structure function in a cone of 10°. However, the difference is still of the same 
order of magnitude as the deviations caused by a variation of 0.1 in a single 
anomalous WW Z coupling. 




-1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 -1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 

cos 9w cos 9w 

Figure 1: The effect on dajd cos 9w of one-loop corrections and anomalous cou- 
plings at -y/i = 174 GeV. Initial state radiation (in a cone of 10°) has been taken 
into account unless otherwise stated. 



7 Summary 

We have presented an event generator for e + e~ — ► W + W~ which includes 

- full 0{<y) corrections, including hard bremsstrahlung, both in the a as in 
the scheme, 

- two-loop leading log effects for the initial state collinear bremsstrahlung, 

- beam polarization (both transverse and longitudinal), 

- polarization of the W's, 
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- anomalous couplings of the triple gauge boson vertices. 

It works in the energy range from threshold to 500 GeV (above this there are 
numerical instabilities, especially in the forward region) . The generation speed is 
of the order of 5 points per second on a workstation. The main limitation is the 
lack of offshell effects; we are working on this problem. 
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A Installation 

The whole package has been written in Fortran 77 with some extensions and has 
been tested on Sur|], DEC CM, NeXT and HP[] workstations. Due to the required 
range in floating point numbers, (m e /mv^) 20 , it does not run as is on non-IEEE 
computers like VAX [] and IBM mainframes. 
The program consists of the following parts: 

- two main programs; one reads its data from a file eeww . dat and calls vegas 
or simplemc to generate the events as a stand-alone program; the other 
demonstrates the use of axmc, which can more easily be tied into standard 
libraries, 

- the switchyard routine wwmc, corresponding to f(x) in Eq. (|l|), which calls 
virt or hard, 

- the main routines with all the physics formulae, 

- the library axo.a, which contains the vegas ||, simplemc, axmc and sup- 
porting routines for integration and event generation, 

- the tensor reduction library aa . a and 

- the library of scalar functions f f . a 21J . 



The whole package can be obtained with anonymous ftp from pss058.psi . ch 
in the subdirectory /pub/eeww.mc, or as a gzip'd tar file. The integration and 
event generation routines can easily be replaced by other packages. 



5 We had problems with the optimizer on Solaris 

6 Do not use +02, change xor to ieor in axo/ranf.f, line 498 and include f lush. hp. o 
7 One could use GFLOAT, but this normally causes problems with libraries 
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B Usage 



There are two ways to generate events. The program eeww reads its input data 
from a file eeww.dat, which it expects in the current directory. This file is read 
a line at a time, with optional lines being skipped over. An example file is shown 
below. 

Standard comment line 

1 order to which the computation is done: or 1 
gmu renormalization scheme: either 'alpha' or 'gmu' 
91.176 mZ Z boson mass 

80.152 mW W boson mass 

130.0 mt top quark mass 

100.0 mH Higgs boson mass 

2 ld-3 IR: 0:normal cutoff, l:Veltman, 2:computed, 3:rescale; cutoff/scale 
2 10. extra initial state coll. radiation: 0:not, 2: included; cone angle 
polarized electron polarization: unpolarized,lefthanded,righthanded, polarized 
0. 0. .9 .9 90. 90. if polarized, polarization matrix (plm,plp,ptm,ptp,f im,f ip) 
polarized W- polarization: unpolarized, polarized, transverse or longitudinal 
polarized W+ polarization: unpolarized, polarized, transverse or longitudinal 
standard anomalous couplings? standard or nonstandard. 

.1.1.1 if nonstandard: dkapg, lamg, xig 

0.1.1.1 if nonstandard: dlz, dkapz, lamz, xiz 
1 number of CMS energies 

200. CMS energies 

simple method: one of 'vegas' or 'simple' 

1 simple: number of points points to generate 

The other method, useful for connecting with other programs, is demonstrated 
by eewwax. This interface consists of an initialisation call, a call which generates 
exactly one event, and an exit routine. The parameters corresponding to the data 
file are passed as arguments. 

The events are output in the routine wweven, which would be the point where 
the W's are decayed and the event analyzed. A sample routine calling jetset 
is included. It has access to arrays of four-momenta, polarization vectors and 
strings, and particle identification codes. Please refer to the comments in this 
files for details. 



C Sample output 

The output of the program with the demo routine wwf eve, given the input file re- 
produced above, is shown below. Please contact the authors if you have problems 
reproducing this or other problems. 
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= EEEEE EEEEE W WW W 
= E E W W W W 

= EEEE EEEE W WW W 
= E E WWWWWW 

= EEEEE EEEEE WW WW 



An order (alpha) generator = 
for W pair production in 
electron positron events. 
J . Fleischer ,F . Jegerlehner = 
K . Kolodzie j , GJ . v . 01denborgh= 



VERSION 1.0 



FF, a package to evaluate one-loop integrals 
written by G. J. van Oldenborgh, NIKHEF-H, Amsterdam 



for the algorithms used see preprint NIKHEF-H 89/17, 
'New Algorithms for One-loop Integrals', by G.J. van 
Oldenborgh and J. A.M. Vermaseren, published in 
Zeitschrift fuer Physik 046(1990)425. 

ffinit: precx = 4.4408920985006D-16 
ffinit: precc = 4.4408920985006D-16 
ffinit: xalogm = 4.9406564584125-324 
ffinit: xclogm = 4.9406564584125-324 
Standard comment line 
eeww: order (alpha) calculation 
eeww: working in the Gmu scheme 
eeww: using masses: 

mZ = 91.176000000000 

mW = 80.152000000000 

mtop = 130.00000000000 

mH = 100.000000000000 
eeww: will compute the cutoff myself 
eeww: including the 2-loop leading log hard and 

exponentiated soft photon effects in cones 

of 10.0000000000000 degrees around the 

beam directions 
eeww: longitudinal polarization electron = 0. 

longitudinal polarization positron = 0. 

transverse polarization electron = 0.90000000000000 

in direction 90.000000000000 degrees 

longitudinal polarization positron = 0.90000000000000 

in direction 90.000000000000 degrees 
eeww: generating all W- polarization states 
eeww: generating all W+ polarization states 
eeww: using standard W couplings 
eeww: running simple event generator 
eeww: CMS energy now is 200.00000000000 GeV 
ranf : using R250 generator 

maxweight increased to 22.968470711021 at event -214 
NP0IN: warning: D4 is not yet supported 
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NPOIN: warning: Bl' seems also not yet supported 
ffxdbd: IR divergent BO', using cutoff 1 . 0000000000000D-24 

ffxdbd: using IR cutoff delta = lam~2 = 1 . 0000000000000D-24 

ffxcOi: infra-red divergent threepoint function, working with a cutoff 

1.0000000000000D-24 
ffzdbd: using IR cutoff delta = lam~2 = 1 . 0000000000000D-24 

maxweight increased to 37.140069173089 at event -212 
maxweight increased to 37.296235986855 at event -4 
simplemc: using as maximum weight: 44.755483184226 

# id E px py pz eps_0 eps_x eps_y eps_z mass 

polarization 

1 11 99.431 0.000 0.000 99.431 0.000 0.000 0.000 0.000 0.001 
POLARIZED 

2 -11 99.431 0.000 0.000 -99.431 0.000 0.000 0.000 0.000 0.001 
POLARIZED 

3 -24 99.117 19.937 2.770 54.723 0.000 -0.138 0.990 0.000 80.152 Y 
TRANSVERSE 

4 24 99.748 -19.937 -2.770 -55.858 0.002 0.931 0.129 -0.342 80.152 X 
TRANSVERSE 

6 22 1.135 0.000 0.000 1.135 0.000 0.000 0.000 0.000 0.000 

7 22 0.001 0.000 0.000 -0.001 0.000 0.000 0.000 0.000 0.000 



The Lund Monte Carlo - JETSET version 7.3 
** Last date of change: 14 Jun 1991 ** 



cross section using weights = 
cross section using or 1 
generated 1 points 
acceptance rate = 100.0000 % 
weight of each accepted event= 
eeww: energy = 200.00000000000 GeV 

c.s. = 18.732903091335 +/- 

chi2 = 0./D0F 



18.732903091335 
44.755483184226 



18.732903091335 



18.732903091335 



total number of errors and warnings 
fferr: no errors 
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